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FOREWORD 


This  report  represents  one  phase  of  a  general  study  of  the 
interaction  of  ions  with  fluid  flow  under  Contract  DA-31-124- 
ARO-D-246,  U.  S.  Army  Research  Office-Durham,  with  Dr.  Henry  R. 
Velkoff  serving  as  principal  investigator. 

The  authors  wish  to  thank  E.  Pejack  for  his  assistance  in 
various  phases  of  the  work  reported  herein. 


ABSTRACT 


Unusually  high  pressure  drops  and  flow  distortions  were 
observed  in  a  previous  experimental  program  involving  laminar  flow 
of  a  gas  in  a  channel  under  the  action  of  a  corona  discharge  in  a 
transverse  electric  field.  A  hypothesis  postulated  by  Velkoff  to 
explain  the  phenomenon  is  extended  to  the  case  of  laminar  bound¬ 
ary  layer  flow  over  a  flat  plate.  The  problem  on  hand  is  found  to 
be  analogous  to  the  laminar  boundary  layer  flow  in  a  transverse 
magnetic  field.  Three  other  mechanisms  proposed  to  interpret  the 
above  experimental  findings  are  also  investigated.  The  increase 
in  viscosity  of  a  gas  because  of  the  ions  is  not  likely  and, 
because  of  the  smallness  of  ion  density,  the  effect  of  ion- 
neutral  particle  interactions  on  the  flow  is  believed  to  be  small. 
One  possible  mechanism  which  may  explain  the  phenomenon  is  the 
secondary  flow  resulting  from  electro-hydrodynamic  instability. 

It  is  found  theoretically  that  Taylor  vortices  can  be  induced 
in  a  quiescent  fluid  between  two  concentric  cylinders  under  the 
action  of  a  corona  discharge.  The  Taylor  Number  of  the  problem 
is  defined  and  shown  to  represent  the  ratio  of  the  destabilizing 
electrostatic  force  to  the  stabilizing  viscous  force.  It  is  also 
found  that  Goertler  vortices  can  occur  in  laminar  boundary  layer 
over  a  flat  plate  provided  the  applied  electric  field  and  the 
charge  density  distribution  satisfy  the  condition  for  instability. 
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INTRODUCTION  AND  SUMMARY 


In  a  previous  exploratory  investigation  into  the  effects  of  ioniza¬ 
tion  on  the  gas  flow  in  a  channel,  interesting  results  were  observed. 
Complete  details  of  that  investigation  may  be  found  by  consulting  refer¬ 
ence  1.  Here  we  shall  only  briefly  describe  these  observations. 

Air  was  passed  through  a  I inch-diameter  pipe  in  which  a  0.004- 
inch-diameter  wire  was  located  concentrically.  A  high  voltage  applied 
to  the  wire  gave  rise  to  corona  discharge  which  provided  ions  in  the  air 
stream.  Under  the  action  of  the  field,  pressure  drops  were  observed  to 
be  doubled,  velocity  profile  distorted,  and  heat  transfer  doubled.  Sim¬ 
ilar  phenomenon  was  also  observed  in  a  5/8-inch  x  5-inch  x  12-foot-long 
rectangular  channel  in  which  ten  parallel  thin  wires  were  located  longi¬ 
tudinally  on  the  center  plane  of  the  channel. 

In  an  attempt  to  explain  the  nature  of  the  phenomenon,  a  hypothesis 
was  put  forth  by  Velkoff  in  Reference  2.  According  to  this  hypothesis, 
an  electric  field  component  in  the  direction  opposite  to  the  flow  is 
induced,  giving  rise  to  a  retarding  electrostatic  body  force.  The  induced 
field  was  found  to  be 


where  u  is  the  flow  velocity  and  K  is  the  ion  mobility  in  the  gas  under 
consideration.  Application  of  this  hypothesis  to  channel  flows  was 
carried  out  in  great  detail  in  Reference  2.  The  close  agreement  of 
theoretical  values  predicted  by  the  hypothesis  with  test  data  aroused 
our  interest  in  extending  the  investigation  to  external  flows. 

The  first  part  of  this  report  describes  an  analysis  extending  the 
above  hypothesis  to  the  case  of  laminar  boundary  layer  over  a  flat  plate. 
An  interesting  analogy  to  magnetohydrodynamics  was  found  for  the  case 
where  the  charge  density  is  constant. 

In  the  remaining  parts  of  this  report,  efforts  are  directed  toward 
the  investigation  of  a  few  mechanisms  proposed  to  account  for  the  phenom¬ 
enon.  The  second  part  presents  a  brief  study  of  the  viscosity  of  a  gas 
mixtrue  consisting  of  ions  and  neutral  particles.  Because  of  insuffi¬ 
cient  knowledge  of  an  exact  expression  for  the  viscosity  uf  ion  gas  at 
room  temperature,  the  discussion  was  based  on  the  data  obtained  for  the 
constituent  gases  at  high  temperature. 

In  Section  III,  the  flow  is  considered  as  a  two-fluid  model.  The 
interaction  between  ions  and  neutral  particles  is  taken  into  account. 

For  the  case  where  the  longitudinal  velocity  component  of  ions  is 
neglected,  the  analysis  resulted  in  an  equation  of  motion  similar  to  the 


1 


one  obtained  by  Velkoff  in  Reference  2.  For  the  case  with  non-vanishing 
longitudinal  velocity  component  of  ions,  the  solution  is  rather  compli¬ 
cated  and  the  effect  of  the  ions  on  the  flow  remains  to  be  determined. 

The  last  part  of  this  report  is  devoted  to  the  study  of  electro- 
hydrodynamic  instability.  It  was  found  that  under  certain  assumptions 
the  quiescent  fluid  between  two  concentric  cylindrical  electrodes  under 
the  action  of  a  transverse  filed  is  governed  by  the  same  eigenvalue 
equation  as  in  the  case  of  couette  flow  between  two  rotating  cylinders. 

A  similar  relation  exists  between  the  case  of  fluid  at  rest  between  two 
infinite  parallel  plates  under  the  action  of  a  transverse  electric  field 
and  that  of  classical  thermal  convection.  It  has  also  been  shown  that 
Goertler  vortices  can  occur  in  the  boundary  layer  over  a  flat  plate  under 
the  action  of  a  transverse  electric  field. 
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I.  LAMINAR  BOUNDARY  LAYER  OVER  A  FLAT  PLATE  IN  A 
TRANSVERSE  ELECTRIC  FIELD 


1.1  Governing  Equations 

In  this  section  we  study  the  effect  of  a  transverse  electric  field 
on  laminar  boundary  layer  over  a  flat  plate.  The  geometricax  configura¬ 
tion  of  the  problem  is  shown  in  Fig.  1.  Let  the  yz-plane  be  taken  as 
the  plane  of  the  boundary  layer  flow  with  z  axis  along  the  plate,  and 
the  y  axis  perpendicular  to  the  plane  wall.  If  p  denotes  the  density 
and  (i  the  viscosity  of  the  fluid,  the  boundary  layer  equation  for  incom¬ 
pressible  flow  in  the  absence  of  an  external  electric  field  is3 


v 


a  d2** 

p  by2 


(1.1) 


Now  suppose  the  fluid  is  positively  charged  through  the  action  of 
corona  discharge  or  the  injection  of  positive  ions  from  external  sources. 
The  charge  density  in  fluid  under  the  suction  of  a  transverse  electric 
field  (Ey),  is  pc  (charge  per  unit  volume).  By  the  hypothesis  postulated 
by  Velkoff,  an  electric  field  (Ez)  opposing  the  flow  is  induced  and  is 
given  by 


Ez  =  -  —  (1.2) 
z  K 

Thus  the  contribution  of  the  transverse  electric  field  to  fluid  flow  is 
in  the  form  of  a  retarding  electrostatic  force.  This  force  is  found  to 
be 


pcEz  =  -  pc 


(1.3) 


Incorporating  this  force  in  the  boundary  layer  equation  (l.l),  we  obtain 
a  modified  boundary  layer  equation  for  the  case  under  consideration 


bw  Pcw 

V  by  bz  pb/2  pk 


(1.4) 


This  equation  together  with  the  continuity  equation 


bf  |  bt 
by  bz 


0 


3 


describe  the  boundary  layer  flow  over  a  flat  plate  under  the  action  of 
a  transverse  electric  field. 


1.2  Constant  Charge  Density 

The  charge  density,  pc,  is  related  to  the  distribution  of  electric 
field  through  the  electrostatic  equation 


V  • 


Pc 

e 


and  is,  in  general,  a  function  of  position.  This  renders  Eq.  (1.4)  dif¬ 
ficult  to  solve.  To  obtain  an  idea  of  what  effects  the  additional  term, 
-Pc/p  w/k,  will  have  on  the  boundary  layer,  let  us  simplify  the  problem 
by  assuming  that  the  charge  density  is  constant.*  Under  this  assumption, 
Eq.  (1.4)  can  be  rewritten  as 


v 


bw  ,  dw  . 

— -  +  w— - +  aw  ■  v  — - 

dz  c*  ~ 


(1.5) 


where  a  =  pc/pK  is  a  constant  and  v  =  p/p  is  the  kinematic  viscosity  of 
tne  fluid. 

1.2.1  Analogy  to  MHD 

The  form  of  Eq.  (1.5)  is  similar  to  that  obtained  for  the  flow 
of  an  electrically  conducting  fluid  over  a  flat  plate  in  the  presence  of 
a  transverse  magnetic  field  fixed  relative  to  the  plate.  The  equation, 
as  was  derived  by  Rossow4  is 


.  dw  .  0U'2  „  2  ^w  /,  *  \ 

v^  +  "S+pHo  w“v^  (1-6) 

where 

0  =  electrical  conductivity  of  the  fluid 
p'  =  permeability  of  the  fluid 


and 


Ho  =  externally  applied  constant  magnetic  field. 

By  comparison,  we  see  that  the  constant,  Pc/K,  corresponds  to  0p’2  Ho2 
in  MHD 


*  Sec. (1.3)  shows  that  this  is  a  good  approximation. 
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1.2.2  Rossow* s  solution 


Equation  (1.6)  was  first  solved  by  Rossow.  He  developed  an 
approximate  solution.  Following  Rossow,  Eq.  (1.5)  is  solved  as  follows. 

Let  us  introduce  the  t  reins  formation 


and  define  the  stream  function  as 

t  83  ^eo  vz  tfo  +  fi  +  mz  f2  +  (mz)3^  f 3  +  (mz)2  f4  +  — ] 

where  woo  is  the  upstream  velocity, 
cc  Pc 

m  =  —  - - and  f0,  fi,  fz,  f3,  f4  —  are  functions  of  11  only. 

Woo  PKWoo 

From  the  definition  of  stream  function,  i.e.,  w  =>  djf/dy  and  v  =  cty/dz, 
v,  w,  ,  d»/bz  and  ZFu/dy2  can  all  be  expressed  in  terms  of  f’s  and 

T).  This  enables  Eq.  (1.5)  to  be  expressed  in  powers  of  mz.  Equating 

the  coefficients  of  equal  powers  of  mz  on  both  sides  of  the  equation,  a 

set  of  ordinary  differential  equations  for  the  f's  are  obtained: 

(a)  2f0"'  =  -  f0"fo 

(b)  2fi'"  -  fo’fi'  -  f0fi"  -  2fxfo" 

(c)  2fa,M  *  2f0'f2'  +  fi'fi’  -  fofs"  -  fifi"  -  3f2f0"  +  2f0' 

(d)  2f3"'  =  3fo'f3'  +  3fa'fi'  -  fofs"  -  2fifa"  -  3f2fi"  -  4faf0”  +  2fx' 

(e)  2f4'"  «  4f0'f4'  +  3fx'f.V  +  2fa'fa'  +  fg'fx'  -  f0f4" 

-  *'.fxf3"  -  3faf2M  -  4f3f1"  -  5f4f0"  +  2fa* 

The  boundary  conditions  v  =  w  =  0  at  y  =  0  and  v  =  dw/3 y  =  0 
at  y  =  00,  when  written  in  terms  of  f's  and  tj,  assume  the  form 

fo  “  fi  *  fa  =*  f3  =  =  —  ■  0  at  i]  =  0 

f0'  =  fi'  =  f2’  =  f3'  =  f4'  —  =  0  at  r)  =  0 

f0'  =1,  f2’  =  -  1  at  t)  =  00 

f  1  '  =  f3'  e  f4'  =  f5'  = - =  0  at  T)  =  00 
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From  the  boundary  conditions  and  Eqs.  (b)  and  (d),  fi  and  f3  can  be  taken 
to  be  zero  throughout  the  flow  field.  Equations  (c)  and  (e)  then  become 


(C)  fa”"  =  f o ' fa '  -  i  f0fa'  -  \  Vo”  +  V » 


(e»)  f4«"  =  2f0,f4‘  +  fa'fa'  -  i  t0Un  -  |  f afa"  -  |  frf©"  +  fa'. 


Equation  (a)  is  the  Blasius  equation  and  Eq.  (c')  was  solved  numerically 
using  Runge-Kutta  method.  The  solutions  were  tabulated  in  Reference  4 
and  are  reproduced  in  Appendix  I. 


1.2.3  Interpretation  of  mz 


It  should  be  noted  that  in  the  above  analysis  mz  plays  the  role 
of  a  controlling  parameter  in  the  flow.  To  investigate  its  physical 
significance,  let  us  non-dimens ionalizo  Eq.  (1.5)  by  introducing  the 
following  dimensionless  variables 


V  -  Re  — 
v* 


TT  V 

W  =  —  , 

Web 


where  L  is  a  characteristic  length  and  Re  -  w«L/v  is  the  Reynolds  number. 
In  terms  of  these  non-dimensional  variables  Eq.  (1.5)  becomes 


c^w  PcL 

V  ~  +  W  —  =  r—K  -  -  W 

dl  cE  dr2  pKwoo 


(1.7) 


Since  all  other  quantities  in  the  above  equation  are  non-dimensional,  the 
quantity  pcL/pKW  must  also  be  non-dimensional  and  can  be  considered  as 
the  product  of  two  non-dimens ional  parameters 


PCL 

pKww 


PcL  .  P 

UK  pWoJj 


We  readily  recognize  that  pw«£/m  is  the  Reynolds  number.  It  is  obvious 
that  pcLa/uK  is  a  new  non-dimensional  parameter  and  is  called  the  charge 
number,  NPc  (Ref.  2).  Now, 


mz  = 


osz 


w„ 


PcZ 

PKwco 


Pcz 

km 


M 

PWaf 


S 


t 


If  we  let  Rez  =  pw^/n  and  NpcZ  =  pcz2/4K,  then  mz  is  the  product  of 
Rez  and  NpcZ.  Further  analysis2  indicates  that  the  charge  number  is 
physically  the  ratio  of  the  electrostatic  force  to  the  viscous  force. 
Meanwhile,  the  Reynolds  number,  Re,  is  known  to  represent  the  ratio  of 
the  inertia  force  to  the  viscous  force.  Thus,  mz  measures  the  ratio  of 
the  electrostatic  force  to  the  inertia  force. 

1.2.4  Velocity  profile 

To  a  first  approximation,  the  velocity  distribution  is  given 
by 


w  =  w#  (f0’  +  mzf2' ) 


(1.8) 


Divided  through  by  w0,  the  free  stream  velocity,  i.e.,  the  velocity  at 
the  edge  of  the  boundary  layer,  on  both  sides  of  the  equation,  Eq.  (1.8) 
becomes 


w 

Wo 


5“  (fo*  +  Bafz1) 


In  Fig.  2,  w/w0  is  plotted  against  q  for  several  values  of  mz. 
The  figure  indicates  that  increase  in  mz  tends  to  retard  the  local  flow. 
From  the  results  of  the  analysis  it  is  predicted  that  further  increase 
in  mz  will  decrease  the  velocity  gradient  at  the  wall,  (dw/dq)__0  • 

When  (dw/dtj)  q  reaches  zero,  corresponding  to  mz  =  0.372,  separation 
will  occur. 


If  the  foregoing  theory  is  applicable  to  this  flow,  it  may  be 
concluded  that  increase  in  electric  field  will  finally  cause  separation. 

1.2.5  Variation  of  boundary  layer  thickness 

If  it  is  desired  to  define  the  boundary  layer  thickness  as 
that  distance  for  which 


7-  -  0-99, 
Wo 


then 


2.  .  Jfl. 

z  Rez 


> 


where  qg  is  the  value  of  q  evaluated  at  y  =  6.  In  Fig.  3,  the  ratio 
6/z  is  plotted  against  Rez.  For  a  given  value  of  z,  the  boundary  layer 
thickness,  5,  increases  with  increase  in  m.  In  other  words,  the  electric 
field  tends  to  thicken  the  boundary  layer.  This  can  also  be  seen  from 
Fig.  2. 
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1.2,6  Local  coefficient  of  skin  friction,  Cf, 


The  local  coefficient  of  skin  friction  is  defined  as 


where  =  |i(6w/<V)s  is  the  shear  stress  on  the  surface.  Using  Eq.  (1.8) 
we  obtain  for  Cf  the  expression 


Cf  =  (0.644  -  I.788  mz  +  ---)  • 

vRez 

This  equation  is  plotted  in  Fig.  4  for  several  values  of  mz.  In  terms 
of  the  charge  number,  Cf  can  be  written  as 


Cf  = 


0.644  _  1.788  N 

n/rcz  (Rez)3^  PcZ 


+ - 


The  variation  of  Cf  with  the  charge  number  is  plotted  in  Fig.  5»  It  is 
seen  that  for  a  given  Rez,  Cf  decreases  linearly  with  increasing  charge 
number.  Therefore,  the  application  of  an  electric  field  reduces  the  skin 

friction. . 


1.3  Variable  Charge  Density 

Consider  that  ions  in  the  fluid  are  provided  by  a  corona  discharge. 
For  the  simple  electrode  configuration  of  two  parallel  plates  (a  plane 
of  fine  wires  is  necessary  to  actually  get  the  corona)  shown  in  Fig.  1, 
the  following  expression  for  charge  density  distribution  is  obtained 
from  Stuetzer's  solution5,  assuming  the  grounded  plate  is  far  away  from 
the  region  of  intense  corona, 


where 

I  =  current, 

A  =  surface  area  of  the  plate, 
E  =  permittivity, 
h  =  distance  between  plates. 
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Assume  that  y/h  <  <  1,  i.e.,  h  >  >  6,  then 


«  1  + 


1  £  + 

2  h  * 


Substituting  this  and  pc  into  Eq.  (1.4)  and  using  Rossow's  method,  we 
finally  obtain 


2f0""  =  -  f0"f0  » 

fa"”  =  fo'fa'  -  £  f0fa"  -  \  faf0"  +  f0'  » 

f4""  =  2f0f4'  +  fafa'  -  i  f0f4"  -  \  fafa"  -  |  f+fo" 


+  f*' +  art  f°  • 

Comparing  these  equations  with  Eqs.  (a),  (o'),  and  (e"),  it  is  found 
that  under  the  assumption  y/h  <  <  1,  only  f4  differs  from  the  case  of 
constant  charge  density.  This  indicates  that  for  this  particular  elec¬ 
trode  configuration  the  charge  density  in  the  boundary  layer  may  be 
assumed  to  be  constant. 
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II.  THE  EFFECT  OF  IONIZATION  ON  VISCOSITY 


2.1  Assumptions 

In  an  attempt  to  determine  the  effect  of  ions  on  the  viscosity  of 
oxygen  gas,  we  consider  a  gas  mixture  of  two  components,  neutral  oxygen 
molecules  and  positive,  singly  ionized  oxygen  ions.  To  facilitate  the 
analys"  ,  the  following  assumptions  are  made: 

(a)  The  mass  of  neutral  oxygen  molecules  is  approxi¬ 
mately  equal  to  that  of  oxygen  ions. 

(b)  Collisions  between  neutral  molecules  as  well  as 
those  between  neutral  molecules  and  ions  are  of 
elastic  rigid  sphere  type. 

(c)  Interactions  between  ions  follow  the  inverse- 
square  law. 

(d)  The  electric  field  intensity  is  small  and  can 
be  neglected. 


2.2  The  Viscosity  for  a  Mixture  of  Neutral  Molecules  and  Ions 

The  analysis  is  based  on  the  general,  expression  for  the  first 
approximation  tc  the  coefficient  of  viscosity  for  a  mixture  of  two  gases 
given  by  Chapman  and  Cowling:6 


[m3i- 


nl2(f 

+H  + 

n-(! 

•**); 

E  .  E  ,  J 

StdiJx  2 [ M2 )x  \ 

nxa| 

'l  +  a'' 

I  ; 

1  t  n21l 

^  +  ikA'1 

O  MX  Rl 

(  E  4A 

1  +  — - : — ; - —  +  ■■ 

- - — t  i  + - +  . . 

[dill  [daJi  2[pi]i[pa]i  3EMxMa 


(2.1) 


where 


Lull 
[uJi 
t  da  Jl 
nN 
n(i> 


=  first  approximation  to  the  viscosity  of  mixture, 

=  first  approximation  to  the  viscosity  of  neutral  oxygen, 

=  first  approximation  to  the  viscosity  of  ions, 

=  number  density  of  neutral  molecules, 

-  number  density  of  ions, 
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nj.2 

mN 

«(i) 

Mi 

Mg 


n 


N 


n2i 


n(i) 


n1 


iN 


mass  of  neutral  molecules, 
mass  of  ions, 

7r— ■  f.\  =  reduced  mass  of  neutral  molecules, 
mN  +  m(i)  * 

m(i) 


mN  m(i) 


=  reduced  mass  of  ions  , 


and  A  and  E  are  constants  determined  by  the  law  of  interaction  between 
molecules.  It  follows  from  assumption  (a),  that 


Equation  (2.1)  then  becomes 


ni2^ 

"Tmli 


E  16a 

2[Ml)l[M2)l  3E 


Under  assumption  (b),  it  is  found  that 


A  =  —  and  E 
5 


-  i  /kml\ >A  1  9 

"2\nj  (dN)  * 


where 


k  ®  Boltzmann  constant, 


aN 


d  =  diameter  of  neutral  molecules  , 


and 


m  =  mass  of  the  molecules  and  T  is  temperature 
Substitution  for  A  from  Eq.  (2.3)  into  Eq.  (2.2)  gives 


(2.3) 


(2.2) 
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U]i  = 


l£>  .16 

I?  012  +  I?  ^ 


E  +  E  +  8_ 
2[mli  2  [  M2  ]  l  15 


16  nig  +  Ito  ngi  +  E 

15  tuili  15  [m.’]i  2[ni]i[Mg]i 


1 6 


-22 
15  E 


(2  JO 


2.3  Evaluation  of  [pi]i  and  [ng]i 

On  the  basis  of  the  elastic  rigid  sphere  model,  the  1  irst 
approximation  to  the  viscosity  of  neutral  oxygen  gas  is  found  to  bec 


[Mx  3j 


(2.5) 


It  should  be  noted  that  all  the  above  expressions  are  obtained  by 
taking  into  consideration  binary  collisions  only.  In  the  evaluation  of 
the  viscosity  of  ions,  difficulty  arises  because  of  the  inverse-square 
law  of  electrostatic  forces.  Since  these  forces  decrease  with  distance 
much  more  slowly  than  the  ordinary  forces  of  interaction,  a  molecule  at 
a  large  distance  from  a  given  molecule  will  also  be  under  the  electro¬ 
static  repulsion  of  many  other  molecules.  Hence,  the  distant  collision 
is  not  binary  but  multiple.  For  the  analysis  based  on  the  binary  colli¬ 
sion  to  be  applicable,  we  can  only  consider  the  special  case  in  which  the 
temperature  is  so  high  and  the  number  density  so  low  that  it  is  suffi¬ 
cient  to  consider  binary  collisions.  This  condition  is  satisfied  at 
6000 °K.  According  to  Chapman  and  Cowling6  the  first  approximation  to  the 
coefficient  of  viscosity,  in  this  case,  is 


[McJi  = 


(^)A- 


(2) 


(2.6) 


with  Ag (2)  given  by 


Ag(2)  =  2  {in  (l+v012)  - 


'J2L. 


1  +  vo1 


where 


and 


d  -  mean  distance  between  neighboring  pairs  of  molecules, 


e  =  electronic  charge  . 
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2.4  A  Numerical  Example 

To  compare  the  values  for  [uiJi  and  [h2  3i>  let  us  evaluate  them  at 
T  =  6000°K.  Assume  the  number  density  ns  of  ions  is  about  1015  cm”3  so 
that  Eq.  (2.6)  applies.  For  this  case, 


and 


=  n  -1/3  a 


10” 5  cm  , 


dkT  _  10”J  x  1.38  x  10“lu  x  6000  „  n 

— p - - rr - -  35*o  » 

e2  2.31  x  10-19 


v0i2  =  (4  x  35. 8)2  =  2.05  x  104  , 

As(2)  =  2  I  In  (1  +  2.05  x  104) - 2-‘P£  x  ±9*  j«17  , 

(  1  +  2.05  x  104) 

[H23i  «  7«1  x  10”6  gm/cm  sec. 

For  reutral  oxygen  molecules  at  T  =  6000°K,  Eq.  (2.5)  gives 


[Hi3i«0.9  x  10”3  gm/cm  sec. 


[Ul3i  and  [te3x  thus  obtained  may  be  substituted  in  Eq.  (2.4)  to 
get  [|i3i  for  given  values  of  n^  and  nli).  However,  without  going  through 
the  calculation,  we  may  draw  the  conclusion  stated  in  the  next  section. 


2.5  Conclusion  and  Discussion 

From  the  above  calculation,  we  see  that  viscosity  of  ions  is  about 
two  orders  of  magnitude  smaller  than  neutral  molecules  at  T  =  6000°K. 

This  indicates  that  the  viscosity  of  the  mixture  is  smaller  than  the 
original  oxygen  gas  of  neutral  molecules.  This  conclusion  is  justified 
by  the  fact  that  in  order  for  the  mixture  to  have  a  maximum  viscosity, 
the  viscosity  of  the  two  pure  components  must  be  nearly  equal.7 

As  was  pointed  out  previously,  Eq.  (2.5)  is  valid  only  when  the 
condition  is  such  that  it  is  sufficient  to  consider  only  binary  colli¬ 
sions.  The  conclusion  made  above  is,  strictly  speaking,  only  valid  for 
this  specific  case.  However,  unless  the  ion  gas  possesses  such  a  peculiar 
property  that  its  viscosity  becomes  greater  than  that  of  its  parent  gas 
at  lower  temperature,  we  may  conclude  from  the  investigation  that  the 
presence  of  ions  intends  to  make  the  viscosity  of  the  mixture  smaller 
than  the  original  gas. 
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presence  of  a  st ro^elec trie °f ield o  fa-6  ^  Gonclusion  even  in  the 
phyrical  reasoning  to  make  the  orpd-in-M  We  d°  nob  ^iave  sufficient 
related  reference!  in  S  S\n°P  ar\We  able  to  find 

problem  deserving  further  investigation^  PP6arS  to  be  “  interesting 
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III.  TWO-FLUID  MODEL 


This  section  describes  an  attempt  to  consider  the  flow  as  consisting 
of  two  interacting  streams,  the  neutral  air  and  the  ion  gas.  For  channel 
flow,  the  neutral  gas  has  only  longitudinal  velocity,  while  the  latter  is 
assumed  to  have  longitudinal  as  well  as  transverse  velocity  components. 

As  we  recall  in  section  I,  the  electrostatic  force  was  considered  to  act 
on  the  fluid  as  a  whole.  It  was  assumed  that  there  was  no  transverse 
velocity  component  in  the  flow.  However,  under  the  action  of  an  exter¬ 
nally  applied  transverse  electric  field,  ions  do  possess  transverse 
velocity.  Thus,  it  appears  logical  to  consider  the  flow  as  the  inter¬ 
action  of  two  streams  of  fluid. 


3.1  Governing  Equations 

Instead  of  including  electrostatic  force  term  on  the  right-hand 
side  of  the  equation  of  motion  for  the  entire  fluid,  we  have  now  to  write 
down  a  separate  set  of  momentum  equations  for  each  species  of  the  fluid, 
taking  into  account  the  change  of  momentum  due  to  collisions  between 
neutral  particles  and  ions.  Such  equations  can  be  derived  from  the 
Boltzmann  equation.  The  detailed  derivation  may  be  found  in  books  on 
plasma  dynamics,  Reference  8,  for  example.  Here  we  give  only  the  final 
expression.  If  the  external  forces  consist  of  an  electric  field  at  id  a 
magnetic  field,  the  momentum  equation  for  the  s^h  charged  species  is 
found  to  be 


(s) 


dt 


+  v»  1 


civ* 


*n(s)  q(s) 


k3  +  n(h»(0  '  *00  (Ei  +  H3k  vJ(")Bk) 


(s) 


=  V 


[vpIJ  -  v-  (s)  ] - i— s  .  ^ 

:nL  i  vi  J  n(s)  °coll 


(3.1) 


where 


Vp  =  itn  velocity  component  of  neutral  particles, 
■p(°)  =  velocity  component  for  the  s*'*1  species, 
n^c)  =  number  density  of  the  s*'*1  species, 

(°)  =  mass  of  particles  of  the  s^1  species, 
pressure-  tensor  for  the  s^1  species 


m 


*..(=) 

Vij 


=  ri(:^m(‘:Kvp(;:)\/  , 
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V^s)  =  ith  peculiar  velocity  component  of  the  s^h  species 
relative  to  its  own  mean  velocity, 

q(s)  =  charge  of  the  s^h  type  particles, 

Ej_  =  i^h  component  of  the  electric  field, 

Gijk  =  the  alternating  unit  vector, 

Bj,.  =  kt  component  of  the  magnetic  field , 

vgn  =  collision  frequency  for  ion-neutral  collisions  , 

Scoll  =  rate  at  which  particles  of  type  s  are  gained  (or  lost) 
because  of  collisions  , 

and  the  <  y  indicates  mean  value  . 

To  apply  Eq.  (3«l)  to  the  problem  under  consideration,  we  make  use 
of  the  following  assumptions: 

(a)  No  external  magnetic  fields  are  applied  and  the  induced 
magnetic  field  is  negligible  because  of  the  small 
current  involved  in  the  problem 

/,  \  The  fluid  consists  of  two  species  only,  ions  and  neutral 
particles.  Collisions  between  ions  and  neutral,  particles 
are  elastic  so  that  SCQ^  vanishes.  Collision  frequency 
Vpjj_  is  constant. 

(c)  The  neutral  gas  flows  between  two  infinite  parallel 
plates  and  the  flow  is  fully  developed. 

(d)  The  flow  is  incompressible 

(e)  The  partial  pressure  of  ions  is  negligible 

(f)  The  ion  gas  is  inviscid. 

(g)  The  positive  electrode  is  the  source  of  ion  gas. 

The  geometrical  configuration  of  the  channel  is  shown  in  Fig.  6. 

If  we  designate  velocity  components  of  the  neutral  gas  in  the  z-  and  the 
y-direction  by  w*'  and  v^.  respectively,  and  the  corresponding  velocity 
components  of  ions  by  w(x)  and  v(0,  because  of  the  assumptions,  vN  =  0 
and,  for  neutral  gas  in  the  z-direction, 

dp 

dx  •  dz  dy2 
J 
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while  in  the  y-direction 


ctyi^  dp 
dx-j  dy 

The  momentum  equations  for  the  neutral  gas  are 

M  =  ^  -  P  vNi(w(1)  -  w1*)  ,  (3.2) 

dy  dz 

=  P  vNiv( i}  ,  (3.3) 

and  those  for  ions  are 

^  -  »„(«"  -  W(i))  ,  (3.4) 

v(i)  ^  E  •  vi*<v(i)  ’  l3-5) 

N  N 

where  p  =  m  n  is  mass  density  of  the  neutral  gas. 

The  nonlinearity  of  Eq.  (3.5)  renders  it  difficult  to  solve  for 

v^1).  However,  the  fact  that  the  transverse  velocity  of  ions  is  the 
drift  velocity  of  ions  in  the  presence  of  an  electric  field,  enables  us 
to  employ  the  concept  of  mobility.  The  drift  velocity  of  ions  is  related 
to  the  electric  field,  E,  by  the  relation 

v^)  =  KE  , 

where  K  is  ion  mobility. 

Equations  (3.4)  and  (3.5)  will  be  solved  for  the  following  two 
cases: 

(a)  Zero  longitudinal  velocity  component  of  ions. 

(b)  Non-vanishing  longitudinal  velocity  component  of  ions. 
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3.2  Zero  Longitudinal  Velocity  Component  of  Ions 

In  order  to  obtain  an  idea  of  how  collisions  between  ions  and  neutral 
particles  will  affect  the  flow,  another  simplifying  assumption  can  be 
made,  namely  negligible  longitudinal  velocity  component  of  ions.  With 
this  assumption  Eq.  (3.2)  reduces  to  the  simple  form 


d^N  _ 
dy2  dz 


pv. 


Ni 


N 

w 


This  equation  is  similar  to  the  one  obtained  by  Velkoff.2 


(3.6) 


The  uncoupling  of  Eq.  (3.2)  from  Eq.  (3.4)  greatly  simplifies  the 
mathematics.  It  is  readily  found  that  the  general  solution  to  Eq.  (3.6) 
is 


w^  =  Ci  cosh^/EZ li  y  +  C2  sinh+fi—^I.  y  -  — i —  —  •  (3.7) 

"  H  *  h  pvfji  dz 

Note  that  in  obtaining  this  solution,  dp/dz  is  considered  constant.  This 
follows  from  the  assumption  that  the  flow  is  fully  developed,  i.e.,  wN  is 
independent  of  z.  The  constants  of  integration  Ci  and  C2  are  to  be  deter¬ 
mined  by  the  boundary  conditions: 

wN(o)  =  w^(h)  =  0. 


The  solution  satisfying  these  boundary  conditions  is  found  to  be 


Integrating  vr  over  the  channel  height',  one  obtains  the  mean  velocity 


Or,  solving  for  the  pressure  gradient  yields 
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3.3  Non- Vanishing  Longitudinal  Velocity  Component  of  Ions 

Now  let  us  examine  the  more  complicated  case  where  the  longitudinal 
velocity  component  of  ions  is  not  negligible.  Equations  (3.2)  and  (3.4) 
must  be  solved  simultaneously.  If  we  let 

1  £P  _  P  P  viN  _  p 

IT  d£  "  C°  ’  —  "  Cl  5  TUT  "  °2  ’ 

where  is  considered  constant  for  simplicity,  then  Eqs.  (3.2)  and 

(3.4)  can  be  rewritten  as 


The  boundary  conditions  are 


wN 


0  at  y  =  h,  w^)  =  w^  =  0  at  y  =  0 


This  set  of  simultaneous  equations  can  be  solved  by  eliminating 
one  of  the  two  variables,  or  by  using  Laplace  transformation.  The 
detailed  calculation  is  given  in  Appendix  II.  Here  we  give  only  the 
final  solution. 


26 


27 


3.4  Determination  of 

From  the  form  of  Eq.  (3.8),  no  explicit  relation  between  the  pres¬ 
sure  gradient  and  the  electric  field  exists.  It  is  through  the  collision 
frequency  for  ion-neutral  particle  collisions,  vjj^,  that  the  electric 
field  manifests  its  influences  on  the  flow. 


Nothing  has  been  said  about  so  far.  The  determination  of 
requires  a  knowledge  of  the  cross  section  for  momentum  transfer  which 
must,  in  general,  be  found  experimentally. 

It  is  shown  in  the  kinetic  theory  of  gases9  that 


where  n(i)  is  determined  by  the  electric  field  and  fjjj  is  independent  of 
the  velocity  and  the  number  density  of  ions  or  neutral  particles.  In 
fact,  fj^j  is  given  by 


e-m*vr2/2kTaj)(vr)vr5dvr  , 
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vr  =  magnitude  of  relative  velocity  between  colliding 
ions  and  neutral  particles 

and 

up  =  transport  cross  section . 

For  rigid-elastic  sphere  model,  fj_N  is  found  to  be 

fiU  -  |  (2m*IIkT)*  diN 

where 

+  <iN) 

and 

d(i)  =  diameter  of  ions, 

N 

d  =  diameter  of  neutral  particles. 

The  relation  between  vNj_  and  is  given  by  m^nN  vNi  =  m^n^Vjjj. 

3.5  Discussion 

No  numerical  data  were  obtained  for  the  analysis  presented  in  this 
part,  for  at  this  stage  of  our  research,  attention  was  called  to  the  con¬ 
sideration  of  stability  problem.  It  is  generally  believed  that,  for  the 
momentum  change  due  to  collisions  between  ions  and  neutral  molecules  to 
have  significant  effect  on  the  flow,  the  number  density  of  ions  must  be 
comparable  to  that  of  neutral  molecules.  Unfortunately,  under  this  condi¬ 
tion  the  assumption  that  the  partial  pressure  of  ions  is  negligible  is 
no  longer  valid. 


IV.  ELECTROHYLRODYNAMIC  INSTABILITY 


As  we  recall  from  classical  hydrodynamic  instability,  flows  between 
two  concentric  rotating  cylinders  will  become  unstable  because  of  cen¬ 
trifugal  forces.  In  the  cases  of  flow  between  concentric  cylindrical 
electrodes  with  corona  discharge,  the  electric  force  in  the  radial  direc¬ 
tion  may  have  a  similar  effect  on  the  flow.  In  fact,  it  can  be  shown* 
that  if  the  components  of  electric  field,  Er,  and  the  charge  density 
gradient,  dpc/dr,  in  the  radial  direction  are  such  that  Er  dpc/dr  <  0, 
an  unstable  situation  could  result.  On  the  other  hand,  Er  dpc/cfcr  >  0 
characterizes  the  stable  equilibrium.  The  fact  that  positive  corona 
between  coaxial  cylinders  does  satisfy  the  condition  for  unstable  equi¬ 
librium**  has  stimulated  our  interest  in  pursuing  investigation  along  the 
line  of  instability. 

The  question  arises  as  to  whether  Taylor  vortices  can  be  created  in 
quiescent  air  under  the  action  of  corona  discharge.  Following  the  clas¬ 
sical  hydrodynamic  instability  problem,  we  propose  to  make  use  of  the 
method  of  small  disturbances. 


4.1  Governing  Equations 

As  before  we  consider  only  incompressible  flow.  The  general  equa¬ 
tions  to  be  used  are  listed  as  follows: 

The  continuity  equation 

V  •  v  =  0  (4.1) 

The  equation  of  motion 

~  +  (v*  •  V)v  =  -  —  Vp  +  vV2  v*  +  —  E*  (4.2) 

ot  p  p 

Equations  from  electrostatics 

(4.3) 

(4.4) 


*  See  Appendix  III. 

**  See  Appendix  IV. 
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Equation  for  current  density 

J*  =  pc(v*  +  KE)  -  DjVpc 
The  conservation  of  charges 

3pc  — > 

+  V*J  =  0 


(4.5) 


(4.6) 


where  v  is  velocity,  Dj_  diffusion  coefficient,  If  current  density  and  the 
other  notations  have  the  same  meanings  as  before. 


4.2  Linearized  Perturbation  Equations 

A  disturbance  of  the  stationary  fluid  will  give  rise  to  the  small 
velocity  (v' )  and  the  pressure  perturbation  (p').  Question  now  arises 
as  to  whether  the  disturbances  will  influence  the  distribution  of  the 
electric  field  intensity. 

Two  approaches  are  possible.  First,  we  may  consider  that  the  per¬ 
turbation  in  the  electric  field  (e^  is  not  negligible.  The  total  electric 
field  is  then  given  by  E  +  e.  The  perturbation  in  charge  density  is 
related  to  e* through  Eqs.  (4.3),  (4.5),  and  (4.6).  It  can  be  imagined 
that  the  electrostatic  force  term  0.1  the  right-hand  side  of  Eq.  (4.2) 
will  give  rise  to  a  rather  complicated  expression  after  being  expressed 
in  terms  of  E  and  5*.  The  formidable  mathematics  thus  introduced  urges 
us  to  seek  an  alternative  approach. 

The  second  approach  to  the  problem  is  to  consider  the  perturbation 
in  charge  density  (pc')  and  assume  the  disturbance  does  not  influence  the 
distribution  of  the  electric  field.  This  assumption  may  be  considered 
to  be  valid  since  the  perturbation  in  charge  density  is  supposed  to  be 
infinitesimally  small  and  the  distortion  of  the  electric  field  caused  by 
it  is  thus  assumed  to  be  negligible. 

With  this  assumption  we  obtain  the  following  linearized  equations 
if  v' ,  p'  and  pc'  are  substituted  in  Eqs.  (4,l),  (4.2),  (4.5),  and  (4.6) 
and  all  quadratic  terms  in  them  are  neglected: 

V  •  v*  =  0  ,  (4.7) 

~  -  Vp'  +  W2  v*  +  -  Pc'E*  ,  (4.8) 

ot  p  p 


=  Pcv>  +  Kpc'E*  -  D^Pc’  ,  (4.9) 


32 


+  v  .  j»  =o  • 
dt 


(U.10) 


Here  only  the  case  of  quiescent  fluid  is  considered.  Detailed  investi¬ 
gations  of  these  equations  will  be  given  in  the  following  sections  for 
two  specific  cases 

(a)  Quiescent  fluid  between  two  concentric  cylinders. 

(b)  Quiescent  fluid  between  two  infinite  parallel  plates. 


4.3  Instability  of  Quiescent  Fluid  between  Two  Concentric 

Cylinders  under  the  Action  of  a  Transverse  Electric  Field 

4.3.1  Formulation  of  eigenvalue  problem 

Let  us  first  consider  the  viscous  fluid  at  rest  between  two 
concentric  cylinders  of  infinite  length  under  the  action  of  a  transverse 
electric  field,  Er  (P).  Let  P,  0,  z  be  cylindrical  polar  coordinates 
and  u' ,  v',  w'  the  corresponding  components  of  the  perturbation  velocity. 
Because  of  the  geometrical,  configuration  under  consideration,  we  assume 
that  the  perturbations  are  rotationary  symmetric,  i.e.,  u',  v' ,  w',  p', 
and  pc'  are  independent  of  0.  Equations  (4.7)  through  (4.10)  for  this 
particular  case  reduce  to 


du'  u'  dw' 

~  +  -2—  =  o 

dr  r  bz. 

9 

(4.H) 

p  Tf  ‘  +  -  fi) 

+  Pc'Er  > 

(4.12) 

9 

(4.13) 

p  =  -  -E-  +  ^w- 

ot  oz 

9 

(4.14) 

-  Di^2Pc '  =  0  , 

(4.25) 

where 

p“  =  mean  charge  density  =  •  E 
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and 


r2  = 


a2 


d 

d? 


d2 

dF 


From  the  form  of  Eq.  (4.13),  we  can  conclude  that  v'  will  be 
damped  out  if  it  is  not  initially  everywhere  zero.  Elimination  of  p' 
between  Eqs,  (4.12)  and  (4.13)  yields 


P 


.  (4.16) 


We  define  non-dimensional  variables 


r* 


r 

ri 


z* 


z 

ri 


T 


TV 

2 

ri 


9 


u* 


9 


Pc* 


Pc'Krg 

V€ 


, 


where 


ri  =  radius  of  inner  cylinder 


and 


€  =  permittivity  . 


Then,  on  dropping  the  stars  on  r,  z,  u,  w,  and  pc,  we  can  write 
Eqs.  (4.1l),  (4.15),  and  (4,l£>)  as 


U  ^ 

T  IS 

=  0  , 

(4.17) 

nKEr 

V 

3Pc  + 

dr 

r?K  _ 
-  pcpc 

_  £i 

V 

^Pc  =  0 

,  (4.18) 

•  (v2  ■ 

-V- 

1  -  r)  • 

ncE, 

:  $££. 

• 

(4.19) 

V 

r  J\bi 

’  dz/ 

pK 

dz 

we  can  make 

the  following 

substitutions: 
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u  =  ui(r)  cos  Az  eaT 
w  =  wi(r)  sin  Az  eaT 
Pc  =  Pci(r)  cos  Az  e0^  > 

where  A  and  a  are  non-dimensional  wave  number  for  the  z -direction  and 
amplification  factor  respectively.  Equations  (4.17),  (4.l8),  and  (4.19) 
then  become 

r  +  ui  +  Ar  wi  =  0  ,  (4.20) 


where  Eq.  (4.21)  is  obtained  by  substituting  for  wi  from  Eq.  (4.20). 

Fc'ir  of  the  six  boundary  conditions  are  that  u  and  w  shall  be 
zero  at  each  of  the  bounding  cylinders.  The  other  two  should  be  on  pc. 
Since  the  perturbation  in  charge  density  is  supposed  to  be  induced  by 
convection,  it  may  be  assumed  to  vanish  at  each  boundary. 

4.3,2  Problem  of  small  spacing 

A  general  discussion  of  the  boundary-value  problem  presented 
by  Eqs.  (4.21)  and  (4.22)  will  clearly  be  a  difficult  matter.  The 
analysis,  however,  simplifies  considerably  when  the  conduction  terms, 
rxKEr/Di  (dpci/dr)  and(r?Kpc/Di€)pcl,  can  be  disregarded  (for  a  justifi¬ 
cation  of  this  approximation,  see  Appendix  V)  and  the  distance  between 
the  cylinders  is  very  small  compared  with  the  mean  of  their  radii,  i.e. 

1  =  r2  -  n  «  \  (r2  =  ri)  =  rQ  .  (4.23) 

We  shall,  therefore,  limit  ourselves  to  this  case.  When  Eq.  (4.23)  holds, 
terms  of  order  of  l/r0  can  be  neglected  and  the  term  l/r  can  be  disregarded 
compared  with  d/dr.  If  it  is  further  assumed  that  the  principle  of 
exchange  of  stabilities  is  valid,  i.e.,  the  situation  in  neutral  stability 
is  given  by  cr  =  0,  we  have,  for  neutral  stability,  the  equations 
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-*( 

askW  , 

,  pK  / 

(4.24) 

= 

dp“ 

— -  Ul 

(4.25) 

eDi 

d? 

With  Er  =£fi(r)  and  dpc/dr  =/gi(r),  where  £and^£  are  con¬ 
stants  having  the  dimensions  of  electric  field  intensity  and  charge  density 
gradient,  respectively,  and,  fi  and  gi  are  dimensionless  functions  assuming 
positive  values  when  evaluated  at  a  given  value  of  their  argument, 

Eqs.  (4.24)  and  (4.25)  can  be  written 


If  we  redefine  non-dimensional  parameters  as 


=  i* 

ri 


Eqs.  (4.26)  can  be  rewritten  as 

(D^-a2)2^  “  aa  —  ^f(S)ocl 

ri  UK 

(D2  -  a2)pcl  =  ~^g(Oui 


where  f  and  g  are  obtained  from  fi  and  gi  by  replacing  r  by  (i/ri  |  +  r0/ri). 
It  is  convenient  to  make  the  transformation 


cDi 


Ul 


Ul 


the  equations  then  takes  the  more*  convenient  forms 
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(D2  -  a2)2Ui  =  -  Ta2fpcl 
(D2  -  a2)pcl  =  gUi  , 


(4.27) 


where 


pDi 


The  boundary  conditions  can  be  written 


Ui  =  DUi  =  PC1  =  0  at  4  =  ±  j. 


Note  DUi  =  0  at  the  boundaries  follows  from  Eq.  (4.20). 

Further  simplification  is  possible,  if  we  consider  the  simple 
case  where  Er  is  constant.  For  this  particular  case,  f  =  1,  £  »  Er, 
g  =  (1  +  Z/r0  0“2  and  A  =  "  €Er/r2.  Because  of  Eq.  (4.23),  g  may  be 
replaced  by  its  value  evaluated  at  |  =  0,  i.e.,  g  =  1.  Eliminating  Ui 
we  obtain 


(D2  -  a2)3pcl  =  -  a2Tpcl>  (4.28) 


The  boundary  conditions  become 


Pci  =  (d2  ”  &2)Pci  =  d(d2  -  a2)Pci  =  0  at  i  =  ±  \  .  (4.29) 


4.3.3  The  variational  principle 

The  characteristic  value  problem  posed  by  Eq.  (4.28)  subject 
to  conditions  Eq.  (4.29)  may  also  be  formulated  using  variational  prin¬ 
ciple  as  follows.10 

Let  P  =  (D2  -  a2)pcl,  and  write  Eq.  (4.28)  as 


(D2  -  a2)2P  =  -  a2Tpcl  . 


(4.30) 


Now  multiply  both  sides  of  Eq.  (4.30)  by  P  and  integrate  over  the  range 
of  | .  The  left-hand  side  of  the  equation  gives 
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r*  */a  ~  Va 

/  P(D2  -  a2)2Pd  i  =  /  PD4Pd  |  -  2a2  /  PD2Pd  i  +  a4  /  P2d  |  . 

*/-Vj  */-'/» 

Since  both  P  and  DP  vanish  at  £  =  ±  \ ,  it  readily  follows  after  two  inte¬ 
grations  by  parts  that 

/  P(D2  -  a2)2Pd  i  =  J  [(D2  -  a2)p]2d  l 
Turning  next  to  the  right-hand  side  of  Eq.  (4.30),  we  have 

f,A  PPoid  5  -  X,  Pc1D2p=id  ?  •  a"  L  Pci2d5  • 

Using  the  boundary  condition  pc  =  0  at  |  =  ±  we  find 


J-'/> 


PpCld  i  = 


f  (DPa  )2d  £  -  a2  f  pc  2d  £  . 

J.v,  -'/»  1 


The  result  of  multiplying  Eq.  (4.30 )  by  P  and  integrating  over  |  is, 
therefore , 


T  = 


f  [D2  -  a2)P)2d  £ 

-'/a 


(4.31) 


^  (DPCl)ad  £  +  a2^  pCl2d|| 
4.3.4  Physical  meaning  of  T 


The  non-dimensional  parameter  T  corresponds  to  the  Taylor 
number  in  flows  between  two  concentric  rotating  cylinders.  It  can  be 
shown  as  follows  that  the  parameter  T  in  the  form 


T  =  - 


uDi 


represents  the  ratio  of  the  destabilizing  electrostatic  force  to  the 
stabilizing  viscous  force. 

Consider  a  perturbation,  of  scale  (40,  in  the  charge  density. 
Associated  with  this  perturbation  is  a  destabilizing  electrostatic  force 
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per  unit  volume  of  scale  (-&<*).  When  convection  takes  place,  convection 
and  diffusion  of  charge  density  have  the  same  order  of  magnitude 
(Eq.  (4.25)),  this  implies  ,hat  the  radial  velocity  component  has  a  scale 
of  Dj/i.  The  viscous  force  opposing  the  convective  motion  is  thus  of 
scale  MDi/^3  per  unit  volume;  this  is  the  stabilizing  force.  The  ratio 
of  the  destabilizing  to  the  stabilizing  force  is 

.  H£g 

MDi 


which  is  T. 


4.3.5  Interpretation  of  the  stability  problem 

The  problem  of  the  solution  to  Eq.  (4.28)  subject  to  the 
boundary  conditions  Eq.  (4.29)  is  an  eigenvalue  problem  leading  to  a 
relation  between  a  and  T,  namely  F(a,  T)  =  0.  It  can  be  foreseen  that 
instability  will  be  possible  if  T  is  greater  than  some  critical  value, 
for  then  the  effect  of  the  destabilizing  electrostatic  force  will  out¬ 
weigh  the  effect  of  the  stabilizing  viscous  force.  The  minimun  critical 
value  of  T  called  Tc,  is  the  minimum  value  of  T,  as  a  function  of  a, 
given  by  F(a,  T)  =  0.  For  T  >  Tc,  it  can  be  inferred  that  a  disturbance 
within  the  band  of  unstable  wavelength  is  an  amplifying,  non-oscillatory 
flow. 


The  variational  formulation  of  the  problem,  Eq.  (4.31),  fac¬ 
ilitates  the  determination  of  an  approximate  value  of  Tc.  If,  for  given 
values  of  a,  T  is  minimized  with  respect  to  variation  of  pCl,  the  resulting 
values  of  T  are  eigenvalues.  In  practice,  a  function  of  pc  is  assumed 
which  contains  one  arbitrary  constant  (-4)  and  satisfies  theHioundary 
conditions,  for  given  a,  the  extremals  of  T,  with  respect  to  variation  of 
are  eigenvalues.  The  minimum  of  these  eigenvalues  for  various  values 
of  a  is  the  desired  critical  value  of  T  (Tc).  The  accepted  value  of  Tc 
is  1707.8  at  a  =  3.13. 

4.3.6  Critical  electric  field  intensity 

From  the  definition  of  the  non-dimensional  parameter  T,  it  is 
clear  that  corresponding  «.o  the  critical  value  of  T  (Tc)  there  is  a 
critical  value  of  the  electric  field  intensity,  Ec,  for  a  given  fluid. 

To  give  an  estimation  of  the  magnitude  of  Ec  for  air,  let  us  consider 
again  the  simple  case  where  Er  is  constant.  For  this  case  the  charge 
density  is  found  to  be 


e 


Er 

T 


and 
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dp- 

15 


eEr 


If  the  gradient  of  the  charge  density  is  evaluated  at  r  =  r0  =  \  (n  +  rz) , 
we  have 


Er  <  0 

dr  tq 


eEv 


and  the  fluid  is  unstable.  Now, 


Tc  =  -  =  1708  - 


MDi 


l-%  rc 


Therefore , 


=  (l708  ^ 


For  air  at  standard  condition 


la. 

i2 


a  = 

Di  = 

e  = 


0.000018  Kg/m  sec, 

0.028  x  10-4m/sec  for  dry  air, 

1.006  x  8.85  x  10”12  ■■g.ga&afr; — 

volt  meter 


If  ro  =  2.5  cm  ana  i  =  0.5  cm  then  r0/ =  10  and  Ec  =  988  volt/m.  Thus, 
for  this  particular  case,  instability  will  set  in  when  the  electric  field 
reaches  the  value  of  9&8  volt/m. 

4.4  Instability  of  a  Layer  of  Fluid  between  Two  Infinite,  Parallel 
Plates  under  the  Action  of  a  Transverse  Electric  Field 

4.4.1  Derivation  of  the  perturbation  equations 

As  a  second  application  of  Eqs.  (4.7)  through  (4.10),  let  us 
investigate  the  case  where  a  layer  of  quiescent  fluid  between  two  infi¬ 
nite  paraxlel  plates  is  under  the  action  of  a  transverse  electric  field, 
E(z).  Let  x,  y,  z  be  rectangular  coordinates  with  z-axis  perpendicular 
to  the  two  plates.  The  plates  are  located  at  z  =  ±  h/2.  To  facilitate 
the  derivation  of  the  perturbation  equation,  we  shall  use  the  notation 
of  Cartesian  tensor  with  usual  summation  convention.10 
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Let  Uj  (j  =  1,  2,  3>  with  1,  2,  3  corresponding  to  x,  y,  z, 
respectively)  denote  the  components  of  the  perturbation  velocity.  Equa¬ 
tions  (4.7)  through  (4.10)  can  then  be  written  as 


cfcxi 

dxi 


=  0 


(4.32) 


dui'  1  dp'  pc  , 

—i  =  -  -  -  +  ~EXi  +  lA^U-! 

ot  p  dxi  p 


(4.33) 


dpc  ,  dps-  i  dE  dpc  , 

—  +  ui  —  +  Kpc  —  ^  +  KE  Xi  t—  -  Di^Pc  «  0  >  (4.34) 

at  ax^  ok±  cbci 

— >  . 

where  A  =  (0,  0,  1)  is  a  unit  vector  in  the  z -direction  and 


V2  »  JL  +  Jf.  +  . 

dxf  dx|  dx§ 


To  eliminate  the  pressure  term  in  Eq.  (4.33) >  apply  the  operator 


curlk  =  eijk 


dX; 


to  the  kth  component  of  the  equation.  Letting 


u-  = 


duk 

€idk-^ 

ox^ 


denote  the  vorticity,  we  have  the  equation 


^i  E  _  dpc  ,  . 

~  -  p  STT ** + 


(4.35) 


Taking  the  curl  of  this  equation  once  again,  we  have 


d  duk  E  .  d2pc 

£1jr  = ; 


dt 


\n  +  vrU- ..  IjS. 

1Jk  “  • 


(4.36) 
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Making  use  of  the  identity 


we  find 


cijk€kAn  3  6i£6jm  “ 


c  ^k 
£lJk  - 


=  e 


ijkekAn 


d^m 


Similarly, 


^p’ 


^Pc 


€ijk€kAn  ^7  ^  =  Ai  SJSI 


-  Ai^Pc  . 


Thus,  Eq.  (4.36)  becomes 


bt  P  \ 


Ai^Pc  -  A 


3^Pi 


dxi^xj 


L) 

'■3/ 


+  vV^ui  • 


Now  multiply  Eqs.  (4.35)  and  (4.37)  by  A^,  we  get 

3S.  =  vV^ 
dt 


and 


.if?!*.  +  f^  +  vVv  , 

at  P  V*2  / 


(4.37) 


(4.38) 


(4.39) 


where  £  =  A^.  and  w  =  A-jU.,  are  the  z-component  of  the  vorticity  and  the 
velocity. 


If  again  we  neglect  the  conduction  terms,  Eq.  (4.34)  becomes 

(4.40) 


a  -  w*  +  Di^Pc 
at  dz 


and  the  required  perturbation  equations  are 


t 


=  vv2^ 
dt 

(4.38) 

8  f  E 

—  V^w'  =  — 

(4.39) 

dt  P 

\5x2  bfz  J 

8pc  _  _  8pc 
dt  dz 

w'  +  D^Pc  . 

(4.40) 

The  fluid  is  considered  to  be  confined  between  two  conducting  plates 
located  at  z  =  ±  h/2.  The  boundary  conditions  on  these  plates  are,  like 
the  previous  case 

£  =  p '  =  w'  =  =  0  at  z  ■■  ±  -  . 

^  dz  2 

4.4.2  Formulation  of  eigenvalue  problem 

Following  the  analysis  in  thermal  convection,  we  now  ascribe 
to  all  quantities  describing  the  perturbation  a  dependence  on  x,  y,  and 
t  of  the  form 


exp[i(kx  x  +  kyy)  +  7  t], 

where  k  =*  >/kx2  +  ky2  is  the  wave  number  and  7  is  a  constant  (which  can 
be  complex).  Suppose  p^,  w' ,  and  £  have  the  forms 


i  =  Ci(z)exp[i(kxx  +  kyy)  +  7  t]  , 
w'  «  w1(z)exp[i(kxx  +  kyy)  +  7  t]  , 
Pc  *  P0j_ )exp[ i(kxx  +  kyy)  +  7  t]  . 
For  functions  with  this  dependence  on  x,  y,  and  t, 


_a 

dt 


=  7 , 


k2,  and  V2  =  -C  -  k2, 
dz2 


and  Eqs.  (4.38),  (4.39)  and  (4.4o)  become 


(4.4l) 


In  terms  of  the  non-dimensional  variables 


(4.42) 


(4.43) 


we  obtain  from  Eqs.  (4.42)  and  (4.43) 

(D2  -  b2)(D2  -  b2  -  a)wL  =  (jjh2^b2  PCi  ,  (4.44) 

if  '  b*  57  °)p=!  =  (57  ^  b2)"1  -  (‘‘•‘•5) 


where  D  =  d/dz*,  b  =  kh  and  cr  =  7h2/v.  The  associated  boundary  conditions 

are 

pc  “0  ,  wi  «=  Dwi  =  0  at  z*  =  ±  |  . 


With  E  =£f(z*)  and  dpc/dz  =^g(z#)  and  making  the  transformation 


we  obtain 


Wi  , 


(D2  -  b2)(D2  -  b2  -  <j)Wx  =  -  b2RfpCl 
(®*  -  62  -  5J  =  BWi  , 


(4.1(6) 


where  R  =  -  is  equivalent  to  the  Rayleigh  number  in  thermal  con¬ 

vection  and  has  the  same  physical  meaning  as  T  in  Section  4.3.  For 
marginal  state  a  =  0,  and  Eqs.  (4.46)  reduce  to 


(D2  -  b2)wi  =  -  b2RfpCi 


(D2  -  b2)pCi  =  gWi 


(4.4?) 


44 


4.5  Instability  of  Laminar  Boundary  Layer  over  a  Flat  Plate  under 
the  Action  of  a  Transverse  Electric  Field 

4.5.1  General  remarks 

It  has  been  shown  in  Section  1.3  that  for  the  parallel  plate 
electrodes  (Fig.  l)  the  charge  density  distribution  in  the  boundary  layer 
is  approximately 


This  indicates  that  dpc/dy  >  0.  But  Ey  <  0,  therefore,  Ey  dpc/dy  <  0 
and  the  flow  is  unstable.  The  situation  here  is  analogous  to  the  cases 
of  laminar  boundary  layers  on  concave  walls  and  on  heated  walls.  For 
the  last  two  cases,  Goertler  demonstrated  that  the  instability  occurs 
in  the  form  of  standing,  longitudinal  cellular  vortices  which  resemble 
the  pattern  of  Taylor  vortices.11}12  It  is  believed  that  this  Goertler 
type  of  vortices  (Fig.  7)  may  also  exist  in  the  problem  under  considera¬ 
tion.  Using  the  assumptions  made  in  the  previous  sections,  we  show  below 
that  this  is  the  case. 

4.5.2  The  equations  of  instability 


Consider  a  flow  past  a  flat  plate  at  zero  incidence;  as  shown 
in  Fig.  1,  the  main  flow  is  directed  along  z,  y  is  normal  to  the  surface, 
and  x  is  along  the  surface  and  perpendicular  to  both  y  and  z;  6  is  the 
boundary  layer  thickness,  and  w,  v,  u  are  the  corresponding  velocities. 


If  all  quantities  are  assumed  independent  of  z,  the  equations 


of  motion  for  three-dimensional  flow  are 


cb  ,  du  . 

¥  +  u5  +  v 


du _ 

p&c 


+  v 


I 


^u 


+ 

*7  ’ 


(4.48) 


k 

at 


+  u 


k  = . 

<v 


(4.49) 

(4.50) 


dx 


0  . 


(4.51) 


And  the  electrodynamic  equations  are 


-  dp" 

Jx  °  pcu  -  Di  — ^  , 

Jy  ■  Pcv  +  KpcEy  -  Di  <^2*  , 

Jz  a  Pcw, 

%  +  u%  +  vT  +  KPc^y  +  KEyly‘  Di(^  +  0  •  (1,<52) 

The  undisturbed  flow  u  =  v  =  0,  w  =  w(y),  p  »  p(y),  Pc  a  Pc(y)  satisfies 
the  above  equations.  Assume,  now,  for  the  disturbed  flow  that 

u  =  ui(y)  sin  a  x  e^ , 

v  *»  vi(y)  cos  a  x  e^  , 

w  =  w(y)  +  wi(y )  cos  a  x  e^  , 

p  =  p(y)  +  Pi(y)  cos  a  x  e^  , 

Pc  88  Pc(y)  +  Pci(y)  cos  a  *  e3t# 


Here  0  is  real  and  denotes  the  amplification  factor,  whereas  X'  =  2K/a 
represents  the  wavelength  of  the  disturbance  at  right  angles  to  the  main 
flow  direction,  ui,  Vi,  wi  are  small  quantities  whose  squares  and  prod¬ 
ucts  can  be  disregarded.  The  vortices  have  the  shape  shown  in  Fig.  7, 
their  axes  being  parallel  to  the  main  flow  direction. 


Assuming  that  perturbations  in  electric  field  intensity  and  in 
electrical  conductivity,  0  ■  KPg,  can  be  disregarded,  we  find  after 
introducing  Eq.  (4.53)  in  Eqs.  (4.48)  through  (4.52) 
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(4.58) 


0Pe  +  Vi 


Eliminating  pi  between  Eqs.  (4.54)  .and  (4.55)  and  substituting  the  value 
of  ui  from  Eq.  (4.57),  we  obtain 

v  -  (8  +  2vo?)  +  c?0  +  va?)vi  -  oP  ££a  Ey,  (4.59) 

dy  djr  p  J 


dSfi 

dy2 


-  (0  +  VQ?)wi  ■  VI  ~ 

dy 


(4.60) 


Di  -  (0  +  Dio?)pc1  -  vi  22S. 

dy2  dy 


(4.61) 


Transforming  the  variables,  let 


i  ■  «(n)  -  a  -  a  b, 


-I-®; 


a2 

-l 


fii, 

Di» 


Wl  -  HqB  Wi, 


n  ■  — 

dq 

where  w©  is  free  stream  velocity.  We  obtain 


(D2  -  c2)(D2  -  c2  -  Olhn  - 


H  Pcx  » 


(4.62) 


(D2  -  c2  -  (TijWi  o  viDW  ,  (4.63) 

(D2  -  c2  -  o-2)pCl  »  vi  ^  .  (4.64) 

Letting  Ey  "  ^(l)>  dp^/dy  B^g(q)  and  making  the  transformation 


6 


Di 


vi- 
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we  obtain,  from  Eqs,  (4.62)  and  (4,64) 


(D2  -  c2)(D2  -  c2  -  ai)Vi  -  -  c2Tfpc1  ) 

(D2  -  c2  -  (TS)pCl  -  gVi  ; 

where  T  *  -  64fijfpvDi  is  the  Taylor  number  for  the  problem, 
conditions  are 


(4.65) 


The  boundary 


Vi  a  DVi,  “  PCl  ■  0  at  tj  ■  0  and  rj  »  »  (4.66) 

Equations  (4.65),  subject  to  (4,66),  define  an  eigenvalue 
problem  for  T,  c,  and  6.  Neutral  stability  is  defined  by  6  *  0. 


4.6  Discussion 

The  theory  of  hydrodynamic  stability  discussed  in  Sections  4.3  and 
4.4  are  primarily  linear  and  valid  for  quiescent  fluid  only.  To  relate 
the  stability  problem  to  the  experimental  findings  cited  previously,  it 
is  necessary  to  take  into  account  the  effects  of  basic  flow  and  finite 
amplitude  disturbance  (the  non-linear  theory).  A  complete  theory  about 
the  interaction  between  the  basic  flow  and  the  finite  amplitude  disturb¬ 
ance  is  as  yet  far  from  being  completely  developed.  However,  certain 
general  understanding  has  been  reached.  Some  of  these  will  be  discussed 
in  this  section. 

According  to  the  linear  theory  a  small  disturbance  will  grow  expo¬ 
nentially  with  time.  As  it  amplifies  it  must  eventually  reach  a  size 
such  that  the  mean  transport  of  momentum  by  the  finite  amplitude  disturb¬ 
ances  is  appreciable  and  such  that  the  associated  mean  stress  (the 
Reynolds  stress 3)has  an  appreciable  effect  on  the  basic  flow.  The 
resulting  distortion  of  the  basic  flow  could  strengthen  the  conversion 
of  energy  from  the  basic  flow  into  the  distrubances.  Since  this  energy 
conversion  is  the  cause  of  the  growth  of  the  disturbance ,  the  rate  of 
growth  of  the  latter  is  filtered.  The  disturbance  is  eLLso  modified  by  the 
generation  of  harmonics  of  the  fundamental  component.  Thus,  there  is  a 
mutual  interaction  between  the  basic  emd  disturbance  peurts  of  the  flow, 
and  a  self -distortion  of  the  disturbfimce. 

When  the  rate  of  conversion  of  energy  from  the  basic  flow  to  the 
disturbance  bdances  the  rate  of  viscous  dissipation  of  kinetic  energy 
by  the  disturbance,  an  equilibrium  state  can  be  achieved  in  which  the 
disturbance  has  a  definite  finite  amplitude  and  the  basic  flow  exhibits 
definite  deviation  from  the  original  laminar  flow.  This  suggests  that, 
in  some  cases,  the  effect  of  instability  is  to  replace  the  origiiml, 
lamineu:  motion  by  another  laminar  motion,  consisting  of  a  basic  motion 
and  a  superimposed  finite  disturbance.  This  may  be  referred  to  as  the 
equilibrium  flow  (or  the  secondary  ilow),  at  a  given  Taylor  number  and 
Reynolds  number,  appropriate  to  the  given  mode  of  disturbeuice .  If  we 
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Impose  the  condition  that  the  rate  of  discharge  for  the  disturbed  motion 
shall  be  the  same  as  the  original,  laminar  motion,  a  larger  pressure  drop 
is  required  because  of  the  work  which  must  be  done  by  the  pressure  gra¬ 
dient  to  maintain  the  disturbance. 


CONCLUDING  REMARKS 


Despite  the  efforts  devoted  to  the  explanation  of  the  phenomenon 
observed  in  the  exploratory  investigation  mentioned  in  the  introduction, 
the  problem  remains  not  completely  solved.  It  appears  that  new  deter¬ 
minable  parameters  must  be  introduced  or  that  one  must  understand  the 
transport  phenomena  of  space  charges  more  completely  than  it  is  known  at 
present.  This  was  borne  out  by  the  fact  that  in  Section  4,  Taylor- 
Goertler  type  vortices  were  shown  to  occur  under  the  assumption  of  negli¬ 
gible  perturbation  in  electrical  conductivity.  The  neglect  of  conduction 
terms  in  the  perturbation  equations  needs  further  justification. 

Lack  of  experimental  data  renders  it  difficult  to  verify  the  analyses 
presented  in  Sections  1  and  4.  Although  no  numerical  calculations  are 
presented  for  the  analysis  in  Section  3>  it  is  believed  that  because  of 
the  smallness  of  the  number  density  of  ions,  the  contribution  of  ion- 
neutral  particle  collisions  to  the  increase  in  pressure  gradient  will 
not  be  important.  The  analysis  of  Section  2  suggests  that  the  pressure 
gradient  increase  is  not  to  be  attributed  to  the  increase  in  the  viscosity 
of  the  fluid  because  of  the  presence  of  ions. 

The  results  of  Section  4  are  of  particular  interest.  They  indicate 
that  theoretically  Taylor -Goertler  type  of  instability  is  possible  under 
certain  circumstances.  It  is  recommended  that  extensive  experiments  be 
conducted  to  verify  the  occurrence  of  Taylor-Goertler  vortices.  If  this 
proves  to  be  successful,  the  next  step  will  be  to  work  out  more  rigorous 
solutions  for  the  eigenvalue  problems  with  various  charge  density  dis¬ 
tribution  and  check  the  critical  Taylor  number  thus  obtained  with  the 
experimentally  measured  values  to  verify  the  validity  of  our  assumptions. 
Because  of  the  highly  mathematical  complexity  and  experimental  difficulty, 
this  is  a  rather  challenging  task. 

Our  analysis  for  flow  between  two  concentric  cylinders  has  been 
restricted  to  the  special  case  of  quiescent  fluid  and  small  spacing. 
Experimentally,  space  charges  are  most  conveniently  provided  in  a  round 
channel  by  a  corona  discharge  from  a  thin  wire  located  concentrically 
serving  as  anode.  To  facilitate  conqjarison  of  theoretical  results  with 
experimental  data,  it  appears  that  the  problem  of  arbitrary  spacings 
should  be  also  considered.  Another  problem  which  concerns  us  is  the 
modification  of  the  analysis  of  Section  4.3  to  the  case  when  an  axial 
flow  is  present. 

In  conclusion,  it  should  be  pointed  out  that  electric  fields  may 
be  used  to  stabilize  laminar  boundary  layers  on  concave  walls  as  well 
as  over  flat  plates,  provided  it  is  possible  to  obtain  a  charge  density 
distribution  such  that  Ey  dpc/dy  satisfies  the  condition  for  stable 
state,  i.e.,  Ey  dpc/dy  >  0, 
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VARIATION  OF  VELOCITY  PROFILES 
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C.  BOUNDARY  LAYER  THICKNESS 
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APPENDIX  II.  APPLICATION  OF  LAPLACE  TRANSFORMATION 


The  simultaneous  equations  to  be 

solved  are 

(H 

+  Cxw^)  =  C0 

(II-D 

-  C2wN  + 

| ^  +  C2^w(l)  =  0  * 

(II-2) 

The  boundary  conditions  are 

w®  =  0  at  y  B  h  ; 

w(i)  =  =  0  at  y  =  0 

• 

The  Laplace  transforms  of  each  term  of  Eqs.  (II -l)  and  (II-2)  are  listed 
as  follows 

Lj^|  -  S»H(S)  -  SwN(o)  -  -S*H(S)-^o  , 

Lfw11}  »  H(S)  , 

L{w(i))  -  H^CS)  , 

»  SH(i)(S)  -  w(i)(o)  -  SH^^S)  . 

After  transformation  Eqs.  (il-l)  and  (II-2)  become 

(S2-Ci)H(s)+  CxH^^(s)  -  %  +  (II-3) 

°  QJy  y»0 

-  C2H(S)  +  (S  +  CajH^S)  =  0  (II-4) 

Elimination  of  H^)(S)  gives,  after  some  algebraic  manipulation  and 
dw1  / 

denoting  - •/  by  B 

dy/  y=0 
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Since  wN  =  0  at  y  =  h,  it  is  found  that 


[(C2eCsh  - 


/2C2Ca  +  Q 
V  Cl 

[(C2ecsh 


/Si  +  cloa  +  OgOahV 
^Ci  C?  Cl  ) 

C2eCsh)  +  (C5eCsh  -  C6eCsh)  -  v/cf  +  4Ci 

Cl 

|ayec5h  _  ec6h)  +  /2a  +  2i|aj(c5eC5h  .  CseCeh) 

-  C2ec6h)  +  (CseCsh  _  c6ec6h)  -  2&  n/c|  +  4Ci] 
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APPENDIX  III.  CRITERION  FOR  THE  ELECTROHYDRODYNAMIC  STABILITY 


Since  the  fluid  is  considered  at  rest,  the  equation  of  motions 
reduces  to  the  simple  form 


PcEr 


d£ 

dr 


for  the  case  of  concentric  cylindrical  electrodes.  This  equation  shows 
that  for  the  fluid  to  he  in  equilibrium  under  the  action  of  an  electric 
field,  the  electrostatic  force  must  be  balanced  by  the  pressure  gradient 
dp/dr.  Consider  now  a  fluid  element  of  charge  density  Pci  at  radius  ri 
and  suppose  that  the  element  is  displaced  to  a  radius  Tz  which  is  greater 
than  ri.  Since  the  charge  density  of  the  fluid  element  remains  unchanged, 
the  new  electrostatic  force  on  it  is  pci  E2,  where  E2  is  the  electric 
field  intensity  at  radius  r2.  The  pressure  gradient  at  radius  r2  supplies 
an  inward  force  of  magnitude  pC2  E2,  Pc2  being  the  original  charge  density 
at  r2.  If  PC2  E2  >  pci  E2,  the  fluid  element  will  be  forced  back  to  its 
original  radius  and  the  fluid  is  said  to  be  in  a  stable  state.  On  the 
other  hand,  if  p  E2  <  pc  E2,  the  fluid  element  will  tend  to  move  farther 
away  from  its  orifinal  position  and  the  fluid  is  in  unstable^  equilibrium. 

Now  if  we  replace  r2  by  r,  rx  by  r-dr,  E2  by  E,  pC2  by  pc,  and  pCi 
by  Pc-(dpc/dr)dr,  where  dr  is  taken  as  a  positive  quantity,  the  condi¬ 
tion  for. stable  state  can  be  written 

ficE  >  (pc  -  ^ 
or 

-  E  dr  <  0 
dr 

Thus,  for  stable  state  E  dpc/dr  >  0.  Similarly,  the  condition  for  unstable 
equilibrium  becomes  E  dpc/dr  <0. 
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APPENDIX  IV.  ELECTRIC  FIELD  AND  SPACE  CHARGE  DISTRIBUTIONS 
FOR  POSITIVE  CORONA  BETWEEN  TWO  CONCENTRIC 
CYLINDERS 


Let  rx  be  the  radius  of  the  inner  cylinder  (or  wire),  r2  the  radius 
of  the  outer  cylinder  (or  tube)  and  J2  the  current  density  at  r  =  r2. 

The  governing  equations  for  the  determination  of  electric  field  and 
charge  density  distributions  are 


dEr  +  Er  _  pc 
dr  r  e 


J  =  PglGSy 


(IV-1) 


(rv-2) 


Combining  Eqs.  (iV-l)  and  (IV-2),  we  obtain 


L  A.  (rE  )  =  zsHsl. 

r  dr  eKrEr 


1  i  i  (r2Ey.2)  =  £S JS.  . 

2  F  dr  '  r  '  eK 


(IV-3) 


The  solution  of  Eq.  (IV-3)  satisfying  the  boundary  condition  Er  =  0  at 
r  =  n  is 


1  _  EL 


From  Eq.  (IV-2)  pc  is  found  to  be 


eJ2r2 

K  p 


(IV-4) 


(IV-5) 


dPc  _  €J2T2 

dr  _~K 


(r2  -  rf): 


Therefore, 


and  the  fluid  is  in  unstable  equilibrium. 


APPENDIX  V.  THE  NEGLECT  OF  CONDUCTION  TERMS  IN  PERTURBATION  EQUATIONS 

The  current  density,  as  given  in  Section  4.1,  is 

J  ■  Pcv  +  KpcE  -  DjVpc  (4.5) 

By  definition,  electrical  conductivity  of  a  fluid  is  given  by 

e  =  p^k+  +  pck“ 

where 


Q  =  electrical  conducitivity 
p*  =  positive  charge  density 
Pc  -  negative  charge  density 
K+  =  positive  ion  mobility 
K  =  negative  ion  mobility 

For  the  problem  at  hand,  pc  =  0,  thus  9  =  p*  K+  and  Eq.  (4.5)  can  be 
written  as 


J  =  pcv  +  0E  -  DjVpc . 

Now,  if  we  assume  that  perturbations  in  electric  field  intensity  as  well 
as  in  electrical  conductivity  are  negligible,  Eq.  (4.9)  assumes  the  form 

^  r) 

J  ■  Pcv  -  DiVpc 


and  Eq. 


(4.15)  becomes 


dpc 

dt 


+  u' 


=  0 


Upon  going  through  the  same  procedure  as  in  Section  4.3,  we  finally 
obtain  Eq.  (4.29). 


67 


_ Unclassified _ 

_ Security  Classification _ 

'  DOCUMENT  CONTROL  DATA  •  RID  ~ 

_ (Steurity  cltnlllctUan  ol  till*,  body  ol  •btlrsci  and  Indotlnf  annotation  uniat  6«  onlotod  whon  tha  wall  ra  port  i*  cl aaal/lad) 

I  ORIGINATIN  C  ACTIVITY  (Corporal*  author)  {a.  NCPORT  SECURITY  C  U  AE»IPIC  A  TION 

The  Ohio  State  University  Research  Foundation  Unclassified 

1314  Kinnear  Road  26  OROUP 

Columbus,  Ohio  43212  _ 


1  REPORT  TITLE 

ANALYTICAL  STUDIES  OF  THE  EFFECTS  OF  IONIZATION  ON  FLUID  FLOWS 


4  DESCRIPTIVE  NOTES  (Typo  ol  topoti  and  inctuaivo  do  fa) 

Technical  Report 


5  AUTHOftfSJ  (Loot  nrnno.  ttrat  n*mo,  Initial) 

Chuang,  T.,  and  Velkoff,  H. 


•  Re  PORT  OATE 

June,  1967 

7#  total  no.  or  PAOEE 

79 

76.  NO  OP  NEPI 

12 

ft*  CONTRACT  ON  GRANT  NO. 

DA-31-124  -ARO-D-246 

6.  PROJECT  NO. 

la  ONIOINATOR’E  REPORT  NUMBERfSJ 

Technical  Report  No.  6 

c 

d 

16.  O 7>l  R^RJ PONT  NOftJ  (Any  otfi.r  man6#rs  Rial  may  6.  aa.IEn.tf 

to  A  V  A  IL  ABILITY /LIMITATION  NOTICES 

Distribution  of  this  document 

is  unlimited. 

II  SUPPLEMENTARY  NOTES 

12  SPONSORINO  MILITARY  ACTIVITY 

U.  S.  Army  Research  Office  -  Durham 

Box  CM,  Duke  Station 

Durham,  North  Carolina  27706 

II  ABSTRACT 

Unusually  high  pressure  drops  and  flow  distortions  were  observed  in  a  previous 
experimental  program  involving  laminar  flow  of  a  gas  in  a  channel  under  the  action 
of  a  corona  discharge  in  a  transverse  electric  field.  A  hypothesis  postulated  by 
Velkoff  to  explain  the  phenomenon  is  extended  to  the  case  of  laminar  boundary 
layer  flow  over  a  flat  plate.  Tne  problem  on  hand  is  found  to  be  analogous  to  the 
laminar  boundary  layer  flow  in  a  transverse  magnetic  field.  Three  other  mechanisms 
proposed  to  interpret  the  above  experimental  findings  are  also  investigated.  The 
ncrease  in  viscosity  of  a  gas  because  of  the  ions  is  not  likely  and,  because  of 
the  smallness  of  ion  density,  the  effect  on  ion-neutral  particle  interactions  on 
the  flow  is  believed  to  be  small.  One  possible  mechanism  which  may  explain  the 
phenomenon  is  the  secondary  flow  resulting  from  electro-hydrodynamic  instability. 

It  is  found  theoretically  that  Taylor  vortices  can  be  induced  in  a  quiescent  fluid 
between  two  concentric  cylinders  under  the  action  of  a  corona  discharge.  The 
Taylor  Number  of  the  problem  is  defined  and  shown  to  represent  the  ratio  of  the 
destabilizing  electrostatic  force  to  the  stabilizing  viscous  force.  In  is  also 
found  that  Goertler  vortices  can  occur  in  laminar  boundary  layer  over  a  flat  plate 
provided  the  applied  electric  field  and  the  charge  density  distribution  satisfy  the 
condition  for  instability. 


5XSI 


Unclassified 

Security  Classification 


LINK  A 

LINK  ■ 

LINK  C 

1  aoi.g 

»  r 

aoua 

WT 

NOLK 

WT 

Electrohydrodynamic  instability 

Taylor  vortex 

Goertler  vortex 

Taylor  number 

Electric  field 

Corona  discharge 

Secondary  Flow 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
of  the  contractor,  aubcontractor,  grantee,  Department  of  De- 
fenae  activity  or  other  organization  (corporate  author)  leeulng 
the  report. 

2a.  REPORT  SECURTY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
"Restricted  Data"  is  included  Meriting  is  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

26.  GROUP:  Automatic  downgrading  is  apecified  in  DoD  Di¬ 
rective  5200. 10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  hsve  been  used  for  Group  3  and  Group  4  as  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  In  all  caaas  should  be  unclassified. 

If  a  meaningful  title  cannot  be  aelected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  In  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  E"ter  the  namefs)  of  authorfe)  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  Initial. 

If  military,  show  rank  end  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE;  Enter  the  date  of  the  report  aa  day, 
month,  year,  or  month,  year-  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7 a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  Le.,  enter  the 
number  of  pages  containing  information. 

76.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  wriu-en. 

86,  8c,  &  8 d.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

96  OTHER  REPORT  NUMBER($):  If  the  report  has  been 
assigned  any  other  report  numbers  ( either  by  the  originator 
or  by  the  sponsor),  also  enter  this  numberfs). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  thla 
report  from  DDC" 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(S)  "All  distribution  of  this  report  is  controlled  Qual¬ 
ified  DDC  uaers  shall  request  through 


If  the  report  has  been  f  rnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

IL  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay 
ing  lor)  the  research  and  development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  s  brief  and  tactual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additional  space  is  required,  a  continuation  sheet  shall' 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS),  (S).  ( C ),  or  (V) 

There  is  no  limitation  cn  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indicstion  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional 


3*0  8SS-SS1 


Unclassified 

Security  Classification 


